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1.  Introduction 


DNA  sequencing  techniques  have  improved  tremendously  in  the  past  few  years. 
An  overarching  objective  of  this  Director’s  Research  Initiative  (DRI)  is  to  explore 
the  use  of  rapidly  improving  state-of-the-art  DNA-analysis  techniques  (commonly 
referred  to  as  “next  generation  sequencing”)  to  begin  to  answer  important 
questions  about  aerosolized  microbes  and  other  biological  particles  in  the 
atmosphere.  Biological  particles  in  the  atmosphere,  such  as  bacteria,  viruses, 
fungal  spores,  and  pollens  (Lighthart  et  al.  1997,  Bowers  et  al.  2011,  Frohlich- 
Nowoisky  et  al.  2012,  Hameed  et  al.  2009),  can  impact  Soldier  health  by  causing 
infections,  asthma,  allergies,  or  other  pulmonary  problems  (Yamamoto  et  al. 
2012).  Additionally,  they  can  impact  weather  by  acting  as  condensation  nuclei 
and  by  absorbing  and  emitting  radiation.  Some  questions  about  the  relative 
abundances  of  airborne  bioparticles  (including  pathogenic  and  allergenic 
organisms)  are  the  following: 

•  How  do  these  vary  with  space,  time,  air  mass  trajectory,  humidity,  and 
temperature? 

•  How  do  they  relate  to  particular  biomes  or  to  local  conditions  (soils, 
vegetation,  and  human  activity)? 

A  long-tenn  objective  is  to  contribute  to  understanding  bioaerosol  threats  to 
Soldier  health  and  develop  monitoring  and  mitigation  strategies  for  bioaerosols. 
Such  strategies  are  needed  to  begin  to  answer  questions  such  as  the  following: 

1)  Are  Soldiers  deployed  to  southwest  (SW)  Asia  exposed  to  significant 
numbers  of  airborne  microorganisms  and  fungal  spores  that  are  different 
from  those  in  the  US? 

2)  How  do  the  abundances  of  the  airborne  microorganisms  and  fungal  spores 
vary  with  location? 

3)  Can  these  differences  between  locations  be  understood,  at  least  partially, 
by  the  climate,  biome,  and  winds? 

Another  objective  is  to  obtain  information  on  the  bioaerosol  backgrounds  in 
different  locations,  which  could  aid  developers  of  biowarfare  agent  detection 
systems.  A  better  understanding  of  how  bioparticles  vary — by  location,  time, 
winds,  etc. — is  needed  for  identifying  strategies  to  minimize  negative  impacts  to 
the  Soldier. 
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The  approach  is  to  do  the  following: 

1)  Collect  aerosol  samples  containing  the  bioparticles  (bioaerosols). 

2)  Obtain  best  estimates  of  the  winds,  air  mass  trajectories,  and 
meteorological  parameters  associated  with  each  sample  using  our  local 
measurements,  local  and  national  databases,  and  the  Hybrid  Single- 
Particle  Lagrangian  Integrated  Trajectory  (HYSPLIT)  and  Weather 
Research  and  Forecasting  (WRF)  modeling  programs. 

3)  Isolate  the  DNA  from  each  sample. 

4)  Selectively  amplify  the  polymerase  chain  reaction  (PCR)  in  the  DNA  near 
evolutionarily  conserved  regions  of  genomes  (the  DNA  coding  for  16S 
ribosomal  RNA  (rRNA)  of  bacteria,  the  18S  rRNA  of  eukaryotes  and  the 
internal  transcribed  spacer  (ITS)  regions  of  fungi). 

5)  Determine  the  sequences  of  these  amplicons  (approximately  6  x  106 
100-base-pair  reads  are  obtained  for  each  of  96  samples  by  collaborators 
at  the  Edgewood  Chemical  Biological  Center  [ECBC]  using  their  Illumina 
Hi-Seq). 

6)  Assign  each  sequence  to  its  most  likely  phylogenetic  categories  (phylum 
to  genus,  depending  upon  the  organism)  using  the  Quantitative  Insights 
Into  Microbial  Ecology  (QIIME)  pipeline  for  taxonomy  assignments 
(Caporaso  et  al.  2012),  cluster  the  sequences,  select  some  of  these  clusters 
to  be  operational  taxonomic  units  (OTUs),  generate  trees  and  networks 
from  the  data,  and  quantify  the  relative  abundance  of  sequences  in 
different  OTUs. 

7)  Use  the  weather  data  and  sequence  data  for  all  the  samples  to  ask 
questions  regarding  how  the  overall  biodiversity  and  the  occurrence  of 
different  sequences  in  the  bioaerosols  varies  and  what,  if  any,  statistically 
significant  variations  are  observed  with  air  mass  trajectories,  source 
location  regions  (e.g.,  Florida  vs.  the  Atlantic  Ocean  vs.  Nebraska),  and 
meteorological  parameters.  An  advantage  of  this  approach  is  that  many 
more  samples  can  be  analyzed  than  with  less  capable  techniques. 

The  data  produced  through  this  effort  could  provide  information  regarding 
differences  in  airborne  microorganisms  in  the  US  vs.  SW  Asia  and  the  Horn  of 
Africa  to  the  US  Army  Public  Health  Command  (USPHC)  and  others  studying 
health  effects  of  deployed  Soldiers.  Understanding  such  correlations  may  help 
lead  to  methods  of  predicting  the  concentrations  of  various  bioaerosols  depending 
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upon  air  mass  trajectories  and  the  meteorological  conditions  along  these 
trajectories. 

The  study  has  only  partially  been  completed.  The  initial  delays  occurred  because 
the  DNA  content  of  the  samples  collected  with  the  Davis  Rotating  Unit  for 
Monitoring  (DRUM)  samplers  was  too  small  to  be  analyzed  effectively  without 
tremendous  amplification.  This  paucity  of  DNA  was  observed  even  with  samples 
collected  over  6-week  periods.  We  were  planning  to  perform  this  study  using  the 
DRUM  samples  in  part  because  US  Army  Research  Laboratory  (ARL)/Battlefield 
Environment  Division  (BED)  has  significant  experience  with  sample  collections 
with  these  DRUMs  in  SW  Asia  and  other  locations.  We  switched  to  collecting 
samples  using  high-volume  (hi-vol)  and  PS-1  samplers  in  order  to  obtain 
sufficient  sample  in  shorter  times.  The  later  delays  occurred  because  of 
difficulties  in  amplifying  ribosomal  DNA  (rDNA)  from  these  complex 
environmental  samples  for  sequencing.  These  issues  have  been  largely  resolved  as 
described  in  Section  4.  Selectively  amplifying  the  16S  rDNA  in  the  presence  of 
excess  18S  rDNA  was  found  out  to  be  more  difficult  than  we  initially  thought. 
That  problem,  combined  with  our  use  of  the  samplers,  which  are  not  size 
selective,  probably  reduces  our  ability  to  detect  bacteria  that  are  in  relatively  low 
concentrations.  We  have  methods  to  further  enhance  the  16S  to  18S  ratio  if  that  is 
desired.  Also,  the  particulate  samplers  routinely  used  by  USPHC  (the  PS-1)  are 
not  size  selective,  and  so  it  is  useful  to  develop  techniques  that  work  with  such 
samples. 

Our  plan  is  to  continue  this  work,  obtain  the  DNA  sequences,  assign  these 
sequences  to  taxa,  and  examine  the  data  for  relations  between  the  taxa  obtained 
and  the  wind  trajectories,  land  or  water  types,  and  atmospheric  conditions. 

2.  Air  Sampling  Techniques,  Samples  Collected  and  Evaluation 
of  Adequacy 

We  first  collected  aerosols  as  stated  in  the  DRI  proposal,  i.e.,  using  the  DRUM 
sampler.  These  samplers  collect  size-  and  time-dependent  samples.  They  have 
been  used  by  BED  and  contractors  to  collect  particles  in  theater.  Samples  were 
collected  using  DRUM  samplers  for  6  weeks  near  Bldg.  202  (Adelphi  Laboratory 
Center  [ALC])  twice  during  the  winter  of  2012-2013.  Several  DNA  isolation 
procedures  were  used  to  obtain  DNA  from  these  samplers.  None  of  the  techniques 
yielded  sufficient  DNA  for  sequencing  reactions,  even  when  all  the  particles 
collected  in  6  weeks  in  2  of  the  8  size  bins  were  combined. 

To  ensure  that  sufficient  DNA  is  collected,  and  is  collected  quickly  enough  that 
variations  in  bioparticle  concentrations  in  times  as  short  as  a  few  hours  can  be 
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examined,  we  collected  particles  using  hi-vol  samplers.  These  samplers  do  not 
differentiate  by  size.  However,  we  found  that  they  sample  fast  enough  to  obtain 
sufficient  DNA  in  times  as  short  as  2  h.  We  collected  124  samples  using  hi-vol 
samplers  between  April  and  February  2014  (Fig.  1).  The  primary  hi-vol  sampler 
was  collocated  with  an  automated  meteorological  sensor  measuring  wind  speed 
and  direction,  temperature,  and  humidity  for  the  majority  of  the  samples  taken. 


Fig.  1  Photographs  of  2  of  the  locations  at  which  samples  were  taken.  Each  shows  the  hi-vol 
sampler  used.  The  upper  photo  also  shows  the  wind,  temperature,  and  humidity  sensor  used. 
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In  February  2014,  we  sent  filters  to  Iraq  to  be  used  by  Soldiers  working  with  the 
USAPHC  to  collect  air  samples  using  PS-1  samplers.  The  filters  reached  the  base, 
but  unfortunately  no  samples  were  taken. 

Because  PS-1  samplers  are  used  by  USAPHC,  we  did  our  last  sets  of  collections 
(approximately  100  samples)  at  ALC  using  these  samplers. 


A  histogram  of  the  sample  length  (in  hours)  for  the  first  series  of  samples  in  2013 
and  February  2014  is  illustrated  in  Fig.  2. 


Sample  Length  (hours) 


Fig.  2  Histogram  of  the  sample  lengths  (h)  used  for  the  first  series  of  samples 


3.  Air  Trajectories,  Meteorological,  and  Air  Quality  Conditions 
during  Sampling  Campaign 


3.1  Air  Mass  Trajectories  Using  HYSPLIT 

Air  mass  trajectories  were  computed  using  HYSPLIT  (Draxler  and  Hess  1997, 
1998,  Draxler  1999).  We  select  this  method  because  it  is  frequently  used  for  dust 
and  smoke  predictions  (Stajner  et  al.  2012),  is  often  used  in  bioaerosol  studies 
(e.g.,  Garcia-Lastra  et  al.  2012,  Sadays  et  al.  2014,  Morris  et  al.  2013,  Urbano  et 
al.  2011  Dannemiller  et  al.  2014),  and  generally  performs  similarly  to  more 
complex  back- trajectory  schemes,  although  more  in-depth  representation  of 
boundary  layer  turbulence  would  be  needed  for  source  area  computations 
(Hegarty  et  al.  2013,  Koracin  et  al.  2011).  Examples  of  air  mass  trajectories 
calculated  using  HYSPLIT  for  2  of  these  samples  (Fig.  3)  illustrate  trajectories 
that  differ  in  source  regions  (southwest  vs.  north  and  east),  time  over  the  ocean 
(none  or  about  1  day),  and  rainfall  (extensive  or  none).  HYSPLIT  trajectories 
were  calculated  for  each  of  the  136  samples  collected  for  the  first  set  of 
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sequencing  runs.  These  trajectories  were  used  as  one  of  the  factors  in  prioritizing 
the  samples  for  sequencing. 

Backward  trajectories  ending  at  2200  UTC  24  Jul  13 
GDAS  Meteorological  Data 
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GDAS  Meteorological  Data 


Fig.  3  Air  mass  trajectories  calculated  using  HYSPLIT  for  2  of  the  samples  collected  with 
the  hi-vol  sampler 
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3.2  Meteorological  and  Air  Quality  Conditions 


3.2.1  Sampling  Campaign  Overview 

Concentrations  of  different  types  of  fluorescent  bioaerosols  are  known  to  fluctuate 
during  the  day.  For  example,  shortly  before  sunrise,  an  increase  in  bioaerosols 
may  be  due  to  nocturnal  spore  emission  (Huffman  et  al.  2012,  Rockett  and 
Kramer  1974).  Others  have  reported  midday  increases  (Hameed  et  al.  2009).  This 
change  may  be  common  to  various  terrestrial  locations  (Lighthart  1997).  As  a 
result,  it  is  important  to  sample  during  different  times  of  day. 

Seasons  also  are  known  to  affect  the  levels  of  bioaerosol  (Jones  and  Cookson 
1983),  and  different  levels  of  biodiversity  as  determined  by  sequences  have  been 
observed  in  different  seasons  (Dannemiller  et  al.  2014).  Meteorological 
conditions  are  also  important,  as  both  the  amount  of  rainfall  (Huffmann  et  al. 
2013,  Schumacher  et  al.  2013)  and  relative  humidity  (Shumacher  et  al.  2013, 
Jones  and  Harrison  2004),  are  correlated  with  changes  in  bioaerosols 
(concentrations  of  different  types).  Thus,  it  is  critical  to  sample  over  a  wide 
variety  of  meteorological  conditions.  The  results  presented  in  the  following 
sections  cover  the  first  set  of  aerosol  sampling  from  April  2013  to  February  2014. 
Meteorological  and  air  quality  data  are  not  available  for  all  samples  during  the 
entire  period,  and  so  the  number  of  actual  data  points  for  each  meteorological  or 
air  quality  parameter  may  vary. 

3.2.2  Meteorological  Conditions 

Relative  humidities  (RHs)  ranged  from  20%  to  above  98%  during  the  sampling 
campaign  (Fig  4).  There  was  no  significant  seasonal  change  in  the  observed  levels 
of  RH,  other  than  that  the  lowest  RH  values  were  observed  during  the  wintertime 
campaigns.  Temperature  ranged  from  above  30  °C  to  -5  °C  during  a  very  cold 
winter  (Fig.  5).  Locally,  winds  were  primarily  out  of  the  west  and  north-northwest 
(Fig.  6).  Winds  from  both  Beltsville  and  College  Park,  Maryland,  were  analyzed 
due  to  uncertainties  in  the  validity  of  wind  data.  Interestingly,  although  local 
winds  were  primarily  out  of  the  west,  synoptically,  air  masses  originated  from 
many  different  locations  and  surface  types  (as  described  in  Section  3.3) 
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Fig.  4  RH  from  the  Adelphi,  Maryland,  area  during  the  sampling  campaign 
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Fig.  5  Temperature  (°C)  from  the  Adelphi,  Maryland,  area  during  the  sampling  campaign 
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Fig.  6  Wind  direction  histogram  from  Beltsville,  Maryland,  (green)  and  College  Park, 
Maryland,  (translucent  gray)  during  the  sampling  campaign 
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3.2.3  Air  Quality  Indicators 

Because  the  majority  of  the  hi-vol  samples  were  collected  from  August  through 
November  2013  and  in  February  2014,  some  atmospheric  contaminants  (ozone, 
PMio)  were  relatively  low  during  this  time  period.  None  of  the  periods  exceeded 
the  Environmental  Protection  Agency’s  (EPA)  24-h  limit  for  either  ozone  or  PMio 
(particulate  matter  with  a  diameter  of  less  than  10  microns).  It  is  relevant  that  the 
highest  levels  of  ozone  are  observed  in  July  through  September,  as  during  this 
period,  both  the  levels  of  oxidants  available  for  reaction  and  the  levels  of  sunlight 
that  cause  the  ozone-producing  photochemical  reactions,  are  the  highest.  In 

-5 

general,  PMio,  as  seen  in  Fig.  7,  was  below  15  pg/m  ,  so  ambient  aerosols,  as 
detennined  by  the  filter  samples,  were  low  from  a  regulatory  perspective. 


Histoaram  of  PM  10 
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Fig.  7  Average  sample  PM10  levels  during  the  sampling  campaign.  The  EPA  National 
Ambient  Air  Quality  Standard  for  PM10  is  150  pg/m3  for  a  24-h  average. 

3.3  Determination  of  High-Resolution  Back-Trajectories  Using 
the  Weather  Research  and  Forecast  Model 

The  influence  of  long-distance  transport  on  the  observed  aerosols  is  uncertain.  To 
better  understand  this  uncertainty,  and  the  potential  factors  which  contribute  to  the 
atmospheric  aerosol  biodiversity,  we  applied  several  modeling  techniques.  We 
investigated  the  calculated  back-trajectories  using  3  different  model  outputs:  from 
a  nested  model  with  a  high-resolution  inner  nest,  a  regional  model,  and  a  global 
model.  The  WRF  (Skamarock  et  al.  2008)  was  integrated  with  nested  grids  with 
27-,  9-,  3-,  and  1-km  horizontal  grid  spacing  centered  over  Adelphi,  Maryland. 
The  9-,  3-,  and  1-km  grids  are  shown  in  Fig.  8.  The  National  Center  for 
Environmental  Prediction’s  (NCEP)  regional  North  American  Mesoscale  model 
(NAM)  data  we  used  have  a  horizontal  grid  spacing  of  12  km.  The  Global  Data 
Assimilation  System  (GDAS)  data  we  used  from  NCEP  have  an  approximate 
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resolution  of  1 1 1x85  km  at  the  latitude  of  Adelphi,  Maryland,  but  is  derived  from 
a  higher-resolution  model. 


Fig.  8  WRF  model  domains  for  the  9-km  (indicated  by  the  extent  of  the  figure),  3-km 
(indicated  by  the  outer  red  box),  and  1-km  (indicated  by  the  inner  red  box)  horizontal  grid 
spacing  nests.  The  27-km  domain  is  much  larger  than  the  3  smaller  domains. 

The  WRF  back-trajectories  were  computed  using  the  Read/Interpolate/  Plot  (RIP) 
program  (http://www2.mmm.ucar.edu/wrf/users/docs/ripug.htnT),  while  NAM  and 
GDAS  back-trajectories  were  computed  using  HYSPLIT.  A  comparison  of  WRF, 
NAM,  and  GDAS  data  for  the  1400  and  2300  universal  time  coordinated  (UTC) 
on  September,  13,  2013,  is  shown  in  Fig.  9.  The  finer-resolution  models  appear  to 
have  more  detailed  features  in  the  back-trajectories.  Although  the  models 
generally  agree  for  the  initial  (near  Adelphi)  trajectory,  there  is  wide  difference 
between  them  during  the  last  few  days  (distant  portion)  of  the  back-trajectories. 
The  end  conclusion  is  that  any  further  work  should  carefully  evaluate  the 
uncertainties,  and  should  use  ensembles  of  back-trajectories  to  ensure  statistical 
significance  of  results. 
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Fig.  9  Partial  back-trajectories  from  the  WRF  (cyan),  NAM  (magenta),  and  GDAS  (black) 
for  1400  and  2400  UTC  on  September  13,  2013,  overlaid  with  land  surface  type 

3.4  Source  Region  Analysis 

During  the  sampling  campaign,  air  masses  at  the  site  had  originated  from  a  wide 
variety  of  locations,  with  many  different  transport  paths.  These  air  masses  passed 
over  different  kinds  of  land  and  ocean  surfaces,  which  may  have  affected  the 
content  and  character  of  the  aerosols  at  the  measurement  location.  To  understand 
the  influence(s)  that  air-mass  trajectory  has,  we  applied  a  time-weighted  residence 
time  analysis  (TWRTA)  using  high-resolution  (250-m)  Geographic  Information 
System  (GIS)  land  surface  data  from  the  North  American  Land  Change 
Monitoring  System  (NALCMS).  Back- trajectories  are  overlaid  with  these  data  for 
illustrative  purposes  (Fig.  10).  To  expedite  analysis,  we  computed  the  percentage 
of  each  number  of  trajectory  sample  points  in  1°  x  1°  bins,  and  used  that  data  to 
compute  TWRTAs,  in  the  manner  of  van  Pinxteren  et  al.  (2010). 

We  discuss  3  samples,  each  with  its  set  of  air-mass  trajectories.  These  samples 
have,  respectively,  the  majority  of  the  influence  from  ocean  (Fig.  10a),  cropland 
(Fig.  10b),  and  forest  (Fig.  10c).  Each  sample  has  approximately  4%-8% 
influence  from  urban  areas.  For  each  of  the  7  trajectories  in  the  “ocean”  sample, 
oceans  contribute  between  36%  and  55%  of  the  TWRTA.  The  surface  type,  which 
contributes  second  most  to  the  TWRTA,  is  either  water  or  the  temperate  broadleaf 
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deciduous  forest,  depending  on  the  individual  trajectory.  During  the  progression 
of  the  trajectory  away  from  the  sampling  location,  initially,  forests  and  cropland 
were  more  important,  but  water  increased  after  a  slight  decrease,  as  did  the  ocean 
influence.  In  summary,  this  sample  is  mostly  ocean,  with  influences  primarily 
from  broadleaf  deciduous  forests,  cropland,  and  coastal  water. 

The  “cropland”  case  has  between  50%  and  60%  cropland  influence,  with 
secondary  influence  (30.8%-42.7%)  from  broadleaf  deciduous  forests.  Other 
minor  influences  (>0.1%)  from  urban,  water,  and  grasslands  are  also  observed. 
The  “forest”  case  is  essentially  an  inverse  of  the  “cropland’  case,  with  54%-55% 
influence  from  broadleaf  deciduous  forest,  39%-39%  cropland,  and  4%  urban, 
with  lesser  influences  from  water,  needleleaf  forests,  and  wetlands. 

Some  general  comments  are  in  order.  Primarily,  back-trajectories  during  these 
samples  spent  the  majority  of  the  time  over  oceans,  cropland,  and  broadleaf 
deciduous  forests,  with  all  other  sources  fairly  minor.  However,  as  urban  areas  are 
located  in  the  immediate  vicinity  of  the  sample  location,  and  the  method  applied 
used  a  linear  weight  (compared  with  an  exponential  weight),  the  influence  of 
locations  which  are  further  away  may  be  inflated. 
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Fig.  10  Superposition  of  back-trajectories  for  three  samples  onto  GIS  data.  High-resolution 
GIS  NALCMS  250-m  data.  Major  contributors  include  coastal  and  inland  water  (dark  blue), 
urban  (red),  cropland  (tan),  wetland  (olive),  temperate/subpolar  broadleaf  deciduous  forest 
(light  green),  and  temperate/subpolar  needleleaf  forest  (dark  green). 

4.  DNA  Extraction  and  Analysis  Techniques  and  Results 


We  developed  a  protocol  for  isolation  of  DNA  from  hi-vol  samples  collected  on 
quartz  fiber  filters.  For  DNA  isolation,  the  filters  were  cut  into  small  pieces  and 
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these  pieces  were  processed  with  a  PowerSoil  DNA  isolation  kit  (MO-BIO). 
Sufficient  DNA  for  downstream  sequencing  reactions  was  obtained  using  this 
approach.  PCR  amplification  was  used  to  determine  if  the  isolated  DNA  was 
amenable  to  standard  molecular  biology  techniques  and  establish  if  16S  DNAs 
and  18S  DNAs  were  present  in  the  samples.  Figure  1 1  shows  an  agarose  gel  used 
to  separate  PCR  amplified  fragments  of  DNA  isolated  from  filters.  Both  the  16S 
DNA  from  prokaryotic  organisms  and  of  the  18S  DNA  from  eukaryotic 
organisms  were  observed  using  18S  primers,  SSU_F04  and  SSU_R22  ,  and  16S 
primers,  357F  and  926R,  as  previously  described  (Sim  et  al.  2012,  Bhadury  et  al. 
2011,  Bik  et  al.  2012).  The  primers  used  for  these  PCR  reactions  amplify  the 
same  region  of  the  16S  and  18S  genes  that  will  be  used  for  organism 
characterization. 
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Fig.  11  Image  of  an  agarose  gel  used  to  separate  PCR  reactions  amplifying  16S  DNA  or  18S 
DNA.  Two  different  preparations  of  DNA  isolated  from  2  different  hi-vol  filters  were  used  as 
templates  for  amplification. 

Ninety-six  samples  were  PCR  amplified  using  the  primer  sets  described  above 
and  shipped  to  ECBC  for  next  generation  sequencing.  Initial  analysis  by  ECBC 
using  an  Agilent  Bioanalyzer  (Fig.  12a)  indicated  the  samples  contained  a  large 
amount  of  DNA  that  was  not  the  appropriate  target  length.  Sequencing  of  these 
samples  would  have  resulted  in  high  levels  of  background,  reducing  the  number 
of  reads  of  the  target  regions.  To  remove  non- target  DNA  from  the  samples  the 
PCRs  were  separated  by  gel  electrophoresis  (Fig.  12c)  and  the  correctly  sized 
bands  were  excised  from  the  gels  and  purified.  DNA  fragments  from  several 
bands  were  sequenced  using  traditional  Sanger  sequencing  to  ensure  the  bands 
contained  the  target  DNA  sequences  (data  not  shown).  Purified  samples  were  sent 
to  ECBC  for  analysis  using  an  Agilent  Bioanalyzer.  The  results  shown  in  Fig.  12b 
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indicated  the  purified  samples  were  enriched  for  the  target  DNA  (length)  and 
were  suitable  for  next  generation  sequencing.  Ninety-six  PCR  amplified-,  gel 
purified- 18S  samples  were  sent  to  ECBC  for  next  generation  sequencing. 


Fig.  12  Analysis  of  the  PCR  samples  for  sequencing  analysis,  a)  Representative  result  from 
the  analysis  of  the  initial  PCR  samples  using  an  Agilent  Bioanalyzer.  The  x-axis  depicts  the 
length  of  the  DNA  and  the  y-axis  is  the  fluorescent  intensity,  which  correlates  to  amount  of 
DNA.  b)  Representative  result  from  the  analysis  of  the  gel-purified  PCR  products  using  an 
Agilent  Bioanalyzer,  c)  Image  of  an  agarose  gel  showing  the  PCR-amplified  bands  for  18S 
RNA  genes.  The  DNA  bands  were  excised  and  purified  for  further  analysis. 


5.  Conclusion 


We  showed  that  to  properly  use  back-trajectories  to  understand  bioaerosol  sources 
ensembles  of  trajectories  with  careful  evaluation  of  uncertainties  were  required. 
We  took  bioaerosol  filter  samples  when  the  measurement  location  was  influenced 
by  a  wide  variety  of  source  locations,  such  as  cropland,  urban,  ocean,  and 
broadleaf  deciduous  forest.  Such  information,  along  with  meteorological  and  air 
quality  data,  will  be  important  in  developing  an  understanding  of  the  factors 
contributing  to  the  bioaerosol  background. 

A  DNA  extraction  technique  was  developed  for  isolating  DNA  from  particles 
collected  by  hi-vol  filters.  Protocols  for  PCR  amplification  of  target  regions  of 
prokaryotic  and  eukaryotic  genomes  were  refined  to  ensure  compatibility  with 
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high  throughput  DNA  sequence  analysis.  After  protocol  refinement,  96  samples 
were  sent  to  ECBC  for  sequencing  on  an  Illumina  Hi-Seq  sequencer  and  we  are 
currently  waiting  for  these  results. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


ALC 

ARL 

BED 

DRI 

DRUM 

ECBC 

EPA 

GDAS 

GIS 

hi-vol 

HYSPLIT 

ITS 

NALCMS 

NAM 

NCEP 

OTUs 

PCR 

QIIME 

RH 

RIP 

rDNA 

rRNA 

SW 

TWRTA 

USPHC 


Adelphi  Laboratory  Center 
US  Army  Research  Laboratory 
Battlefield  Environment  Division 
Director’s  Research  Initiative 
Davis  Rotating  Unit  for  Monitoring 
Edgewood  Chemical  Biological  Center 
Environmental  Protection  Agency 
Global  Data  Assimilation  System 
Geographic  Information  System 
high-volume 

Hybrid  Single-Particle  Lagrangian  Integrated  Trajectory 

internal  transcribed  spacer 

North  American  Land  Change  Monitoring  System 

North  American  Mesoscale 

National  Center  for  Environmental  Prediction 

operational  taxonomic  units 

polymerase  chain  reaction 

Quantitative  Insights  Into  Microbial  Ecology 

relative  humidity 

Read/Interpolate/  Plot 

ribosomal  DNA 

ribosomal  RNA 

southwest 

time-weighted  residence  time  analysis 
US  Army  Public  Health  Command 
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UTC 


universal  time  coordinated 


WRF 


Weather  Research  and  Forecasting 
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